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\Q ' We have studied the effect of anisotropies on the quantum phase transition of the Kondo necklace 

' model in dimensions D=l, 2 and 3. Both the anisotropy 8 of the inter-site interaction term and 

anisotropy A of the on-site Kondo interaction have been included. We use a bond operator method 
with constraints implemented in mean field approximation. Starting from the paramagnetic phase 
we determine the critical ratio (t / J) c of the quantum critical point and associated scaling exponents 
of the Kondo-singlet gap. We show that in the case of easy-axis type anisotropy S > 1 a qualitatively 
new behavior in comparison to the conventional Kondo necklace model with (<5,A)=(0,1) appears. 
We have also obtained the antiferromagnetic order parameter in the long range ordered phase for 

t > t c . 



o 



PACS numbers: 75.10 Jm, 75.30 Mb 
£j ■ I. INTRODUCTION 



I 



OS 

in 



Quantum phase transitions in heavy fermion systems between a Fermi liquid and antiferromagnetically (AF) ordered 
state have been at the focus of research recently 0, This is due to mainly two reasons: Firstly in the vicinity of 
the quantum critical point (QCP) which is usually reached by applying external or doping-induced chemical pressure, 
the intrinsic energy scale of the compound, i.e., effective Fermi temperature T* or Neel temperature Tjv vanishes. 
This leads to a breakdown of the Fermi liquid picture and non-Fermi liquid behavior of thermodynamic and transport 
quantities is observed close to the QCP Q. It is therefore of great interest to understand the approach to the 
i quantum critical region within suitable theoretical models. The most important among them is the Kondo lattice 

model consisting of a free conduction band and an on-site AF Kondo interaction which favors nonmagnetic singlet 
CN| . formation. In the second order it also leads to the effective inter-site RKKY interactions which favor magnetic order. 
Their competition leads to the appearance of the QCP. In such a picture only spin degrees of freedom are involved 
in the quantum phase transition. Therefore the Kondo lattice model may be replaced by a simpler model where the 
itinerant hopping part is simulated by an inter-site interaction of the itinerant spins. In the work of Doniach [J] it was 
shown that in dimension D=l this replacement, leading to the 'Kondo-necklace' model is indeed exact assuming that 
the inter-site interaction is of the XY-type. Later more general anisotropic inter-site interactions have been considered 
\Q ' including the z-component of conduction-electron spins |5| . In the itinerant picture this would be equivalent to adding 
an interaction term for conduction electrons, i.e., considering Kondo spins that are screened by correlated conduction 
electrons 0. In addition the Kondo-necklace model was considered in higher dimension 0. In this case the direct 
connection to the original itinerant Kondo lattice model is lost. One then has to consider the Kondo necklace model 
for D > 2 as a model in its own right. In fact exact diagonalisation results j9j for finite clusters show that the two 
models are still closely related even for D = 2. 

In the present work we want to study the possible quantum phase transition in the D = 1-3 Kondo-necklace type 
model under rather general assumption of the anisotropy of both the inter-site interaction and on-site Kondo terms. 
£J . They are characterised by the pair of parameters (8, A). For the former 8 ^ describes the deviation from the XY- 
type interaction associated with free conduction electrons in ID and for the latter A ^ 1 describes deviation from 
the isotropic on-site Kondo interaction. The A anisotropy is always present in real Kondo compounds like Ce-based 
■ intermetallics due to the crystalline electric field (CEF). 

For uniaxial symmetry it splits the J = | multiplet into a sequence of Kramers doublets. The ground state doublet 
|±) may be described by a pseudo-spin S=i. Projecting the Kondo term to this subspace then leads to different on- 
site exchange for S XjV and S z components. This defines a local anisotropy ratio A = 3 Z /3 X of the pseudo-spin Kondo 
term. In real Kondo compounds this ratio may assume any value between the easy axis (Ising-like) case A — > oo 
and the ea sy-p lane (xy-like) case A = 0. Examples of tetragonal compounds are CeRu2Si2 for the former and 
YbRh 2 Si 2 BjJ for the latter. 

We study the quantum phase transition from the paramagnetic (Kondo-singlet) side as function of the control 
paramter t/J giving the ratio of inter-site to on-site interaction strength and as function of the ansisotropies (5, A). 
Previously this has only been performed in the special case (<5,A)=(0,1) 0. For the case of general (8, A) we use 
the same bond operator representations of spin variables where on-site constraints are implemented in mean field 
approximation. 
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FIG. 1: The energy gap E g versus hopping strength in the one dimensional lattice. Different plots show the effect of local 
anisotropy, < A < 1. 



In Sec. [H] we give a brief definition of the anisotropic Kondo necklace model and the bond-operator transformation. 
In Sec. IIIII the ground state energy and excitation spectrum is calculated in mean field approximation. Using the 
numerical results the quantum phase transitions and phase diagrams are discussed in Sec. II V I The extension of the 
mean field approach for the antiferromagnetic ordered phase is presented in Sec.0 Finally Sec. IVII gives the summary. 



II. ANISOTROPIC KONDO NECKLACE MODEL 



We investigate the anisotropic Kondo necklace model which is defined by the following Hamiltonian 

H = H t +Hj=tJ2 (TfTf+TfTj+STfrf + J^TfSf+TfSl+ATfSf). (1) 
<*iJ> * 

In the above Hamiltonian, r" represent the a-component of spin of the itinerant electron at site n and S" is the 
a-component of localized spins at position n. For the exchange coupling between the itinerant and localized spins we 
have used J x = J = 1 as the reference energy scale in all figures. The hopping parameter of the itinerant electrons 
is proportional to t with the anisotropy in the z-direction given by 5. Therefore the above model has three control 
parameters: t/J and the anisotropies {6, A). There are two main reasons for considering the anisotropy both in 
the itinerant and localized parts of the interactions: (i) In real materials where the Kondo effect has been observed 
there exists the anisotropy in the interaction between the localized and itinerant electrons, (ii) The anisotropy in the 
itinerant part of the interaction enables us to investigate the effect of symmetry on the results of the Kondo necklace 
model compared with the Kondo model. In the present model the case of 8 ^ 1 has U(l) symmetry and for 6=1 
it has SU(2) symmetry, the same as in the genuine Kondo model. Moreover, in one dimensional lattice nonzero <5 
represents the interaction between itinerant electrons [l2| . 

Our approach is based on the strong coupling limit J/t — > oo(t = 0) where the model is composed of independent 
pairs of spins (r, S). The Hilbert space of a pair of spins consists of 4 states where it can be represented by a singlet 
and a triplet. This basis can be created out of the vacuum by singlet and triplet creation operators: \s) — s^O) 
and \t) — ijj|0) (a ~ x,y,z). In terms of singlet-triplet operators the spin operator of the localized and conduction 



3 




FIG. 2: The energy gap E g versus hopping strength in the two-dimensional lattice for 5 = 0. The inset shows the log-log plot 
of Eg versus \t — t c \ for A = 0, 1 where the gap scales like E g ~ \t — t c \ v with v ~ 1 close to the critical point. 



electrons are given by 

1~n,a — ^( ^ri^n^a ^n.a^n ^af3^^ n ^^n.'y ) ? (2) 

where (a, /3, 7) represent the (x, y, z) components and e is the totally antisymmetric tensor. This type of representation 
has been introduced first by Sachdev and Bhatt |l3j called bond operators. The bond operators satisfy bosonic 
commutation relations [s n , s|J = 1, \t n ,ait n r\ = ^a,0 and [s n ,t„J = 0. The physical states are obtained by the 
local constraint sJjS„ + ^ Q in.a = 1- Using Eq.© the Hamiltonian (Eq.Q) may be expressed in terms of bond 
operators. As discussed in the next section it then reduces to a simple form within the mean-field approximation. 



III. SINGLET CONDENSATION AND GROUND STATE PROPERTIES IN MEAN-FIELD 

APPROXIMATION 

We have followed the mean-field approach introduced in Ref0 where the local constraint is replaced by a global one 
for the average amplitudes (s n ) and (t n _ a ). In this respect we start from the strong coupling limit J/t -^00 (t = 0) 
where the ground state is composed of the direct product of local singlets. In this limit we have (s„) = (s^) = 1 and 
(t n ,a) — (it, a ) ~ where the ground state represents the pure condensation of local singlets. Turning on the hopping 
of conduction electrons (t ^ 0) the triplet occupation at each site becomes nonzero but still very small as will be 
shown later. Thus, for a small value of t/J we consider a mean field value for (s„) = (sjj) = s and the excitations 
above this background are defined by the triplet states which define the mean field Hamiltonian. The hopping term 
of the Hamiltonian which is modeled by the XXZ interaction between the spin of conduction electrons will be written 
as a sum of three terms, i.e., Ht = H 1 + H2 + H3 associated with the different types of singlet and triplet interactions. 
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FIG. 3: The occupation of singlets in the ground state (< s >= s) versus the hopping strength up to the critical point t c for 
two different values of anisotropy A = 0, 1 in the two dimensional lattice. The crosses (x) indicate the critical point for each 
case. The deviation from fully occupied state (< s >= 1) is very small for t <t c . 



The first term is given by 

Hi = J Y ( S i t nA S la t m^+tl lia S m )+h.c}j (3) 

<n,m> a.—x,y 

+ J H (4 l i™, Z ( s L i m,2+4, Z Sm) + ^- c -)- (4) 
<n,m> 

The explicit form of H2 and H3 are given in the appendix. However, H2 contains the interaction between the triplet 
bosons which will have negligible effect on the result of considering only Hi because the average triplet occupation is 
small. The contribution of H 3 to the ground state energy is identically zero in mean field approximation, since each 
term of H 3 consists of three triplets and one singlet operator. 

The exchange term (Hj) is diagonal in terms of the bond operators, 

Hj = j J2 ( - ( 2 + A ) s « s ™ + ( 2 - A )C*«.* + A(tJ, >x t„, x + . (5) 

n 

The physical constraint is imposed by adding a Lagrange term at each site with an associated chemical potential /U„ 
to the Hamiltonian 

H = H 1 +H J + ^2n n (s$ l s n + ^2 t n , a tj,,a-l)- ( 6 ) 

n a—x : y : z 

The mean field Hamiltonian, H m f — (H) is obtained by taking (s n ) = (s+ ) = s and a global value for the chem- 
ical potential, /U„ = fj,. After performing the Fourier transform and the Bogoliubov transformation the mean field 
Hamiltonian is written in the following form, 

H mf = E +^2 u a{k)a\ a a k , a . (7) 

k a=x,y,z 

In this equation, Eq is the ground state energy, u> a (k) is the excitation energy of the new bosons defined by the 
number operator a\ a ak, a - The new bosons can be expressed in terms of the bond operators via the Bogoliubov 
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transformations, 



where k is the momentum and 



a k ^ a = cosh(9 k , a )tk, a + sinh(6'fe !Q )tL fe!Q , 
a'_ ka = smh(0fc, a )i fc , Q + cosh(6» fc!Q )tL fc a , (8) 



9 _2/.(fc). _2^(fc) 



The functions f a (k) and d a (k) are defined in terms of the coupling constants by 



/»(*) = fy(k) = ^7(fc); dx(fc) = d y (k) = M + f + ^-7(fc), 



/«(*) = ^(k); d z (k) 2 Jj—l± + ^(k), (10) 

where J z — J A, j(k) = Yli=i cos(fci) and D is the dimension of rectangular Brillouin zone (BZ) associated with the 
D-dimensional rectangular lattice. The bosonic excitation energies in Eq. J7J) are obtained as 



^W = V#F1M, a = x,y,z. (11) 
The ground state energy in D dimension is then given by 



E« = n(^-1)-{^±±)?)+\Y. E (««(*)-*.(*))■ (12) 
Minimization with respect to /i and s implies 



k a=x,y,z 



^-0 ^-0 (13) 



Explicitly this leads to 



1 x - djjt) = 5 _ _ 2 
2N ^ uJk) K 2 ' 

koc 

^ Li — ^ — T«( fe )] - 2J ^ 1 + o~tJ - j-J- ( 14 ) 

ka 



Here we defined J x ,y(k) — 7(A) and 7z(fc) = 5~f(k). These equations will be solved selfconsistently for the average 
Kondo singlet amplitude s and the chemical potential fi. It is worthwile to write the m.f. Hamiltonian Eq.Q in a 
different form as 



Hmf = Eq+Y^ E W «(fc)(4,a a fe,«+ 9) 



2' 

k a—x,y,z 
' *2J X -f- J z , _2 1^ ,5 _2> 



E = -N[^^{s' + -)+^-s')) (15) 

Because s < 1 the first part (£? ) is always negative and represents the m.f. singlet condensation energy. The second 
part is the (positive) energy of excited triplet bosons. The ground state has only a contribution from zero point 
motion. Both m.f. condensation energy and energy of zero point quantum fluctuations contribute to the total energy. 
Their balance determines the critical t c for the quantum phase transition where the zero point fluctuation amplitude 
of bosons diverges, i.e., uo a (Q) becomes soft. 



IV. NUMERICAL RESULTS: THE QUANTUM CRITICAL PHASE DIAGRAMS 

We will now obtain the numerical solutions to the coupled equations (|13fl . For this purpose we first consider the 
XY- case, i.e., 5 = and arbitrary anisotropy A for the Kondo term in Eq. Q). We remind that in D=l the XY-case 
is equivalent to the genuine Kondo lattice model with free itinerant electrons of a band- width W=4t. Later we will 
also investigate the case of nonzero 5. 
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FIG. 4: The energy gap E g versus hopping strength in the three dimensional lattice. The gap scales like E g ~ \t — t c \ v with 
^~0.55. 



A. XY-case: <5 = 

In this case the z-polarised branch of excitations has a dispersionless value u> z (k) = luq. However, the other two 
branches of excitations (lo x = u> y ) show a dispersion which has a minimum at the antiferromagnetic reciprocal vector 
at the corner of the BZ, e.g., Q=(7r,7r) for D=2. The minimum value of excitations defines the energy gap 

£ »= ( ''+T>/^lbi>- (16) 

In the present case the energy gap depends on the two parameters t and J z . The gap E g defines the energy scale 
for the Kondo singlet phase. When this energy scale approaches zero the other type of excitations become important 
and the model will encounter a phase transition from Kondo singlet to the antiferromagnetic phase. We will see that 
there exists a critical value for t = t c where this phase transition occurs in D = 2, 3 while no transition happens in 
D = 1. This is similar to the results in Ref. [D for the special case (5,A) = (0,1) and to the previous results of MC 
simulations 8] for the ID model. We have solved the mean field equations l|13|) numerically for general values of the 
coupling constants (i, J) and anisotropies (<5, A). The resulting value of \i and s are replaced in Ea. (|16|) to obtain the 
energy gap and track the location t c where it vanishes. 

In Fig^we have plotted the energy gap versus t (the hopping strength) for different anisotropies < A = if- < 1 
in D = 1. The value of gap is reduced by decreasing the anisotropy from the isotropic (A = 1) case to the XY-type 
(A = 0) Kondo interaction. However, it is always nonzero and never vanishes. Thus, the effect of anisotropy does 
not change the universality behavior of the D — 1 case where the system is always in the Kondo singlet phase. We 
have also considered values with A > 1 (easy axis- or Ising-type regime) and found a similar behavior as in Fig^ 
provided we change the energy scale from J x to J z . It can be seen from the expression for the energy gap that there 
is a duality between the < A < 1 and A > 1 regions by interchanging J x < — ► J z . For theoretical reasons we might 
also consider the negative values of anisotropy (— 1 < A < 0). Its behavior is similar to Fig^ while the minimum 
of gap is reduced and vanishes at t = for A = — 1. Actually, at A = — 1 the energy difference of the local singlet 
and the two triplet states becomes zero which leads to the vanishing of the energy gap. However, in this limit the 
above mean field approach is not reliable where the triplet contribution can not be neglected in the ground state. For 
negative values of A the average value of the triplets is comparable with the singlets and should be taken into account 
on the same footing. 
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FIG. 5: The energy gap versus t for different values of S in the one dimensional lattice. The gap is always nonzero for 8 > 0. 



Let us now present the results for the square lattice (D = 2) shown in Fig|2| We have plotted the energy gap versus 
t for different values of < A < 1. For all values of A the gap vanishes at the critical point t c , where the transition 
from Kondo singlet to the antiferromagnetic phase occurs. The dependence of the critical point on the anisotropy is 
linear and is fitted well by t c ~ 0.35J(A + 1) which is also valid for negative values of — 1 < A < but has not been 
shown here. We will derive an equation for the critical temperature of D = 2, 3 in terms of the coupling parameters 
in Sec. The linear dependence of t c on A shows that the anisotropy term is important to stabilize the Kondo 
singlets. Thus, for higher A, higher value of the hopping term is required to destroy the Kondo singlet and induce the 
antiferromagnetic order. However, the qualitative behavior of the energy gap versus t is the same for all A. Our fine 
tuned data plotted in the inset of Fig|21show that the gap vanishes close to t c with a critical exponent v according to 

E B ~\t-t D \ v , (17) 

where v ~ 1 for D = 2. The exponent is the same for all anisotropics and this is in agreement with the observation 
that the local interaction terms will not change the universal behavior. The log-log plot of E g versus \t — t c \ for all A 
fall on each other which justifies the scaling behaviour close to critical point. The extension of our calculation to A > 1 
(Ising regime) shows qualitatively similar results. We have also plotted the value of s —< s > versus the hopping 
term in Fig. [3] for two different values of anisotropy, A = 0, 1 in the case of two dimensional square lattice. These 
data justify our assumption that the expectation value of singlets in the ground state is close to one up to the crit- 
ical t c value. It supports our claim to neglect the occupation of triplets in the ground state in the mean field approach. 

Our results for D — 3 are shown in Fig. The gap versus t shows similar behavior qualitatively for all values of 
A. There exists always a critical point (t c ) which shows the transition from Kondo singlet to the antiferromagnetic 
phase. The critical point depends on the anisotropy linearly, t c ~ 0.188J(A + 1). The exact expression for the critical 
point is given by an equation in Sec[V] However, the dependence is weaker than in the case of D = 2. The reason is 
related to the increase of the number of nearest neighbor sites which allows hopping more easily. Thus, the hopping 
term is more effective here. The analysis of log- log plot of Fig. 0] close to the critical point gives the gap exponent to 
be v ~ 0.55 for the three dimensional cubic lattice. 
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FIG. 6: Two dimensional lattice: (a) The behavior of gap for two different regimes of 5 versus t in two dimensional lattice for 
A = 1. The scaling of gap close to the critical hopping (t c ) is obviously different for < S < 1 and S > 1, (S — 2,3,5). (b) 
Staggered magnetization (m s ) versus t in the antiferromagnetic phase, just above the quantum critical point at t c . The local 
anisotropy is A = 1 and the itinerant one varies within < 5 < 1. 



B. general ansisotropic case 5^0 

Our results in the previous case, S = 0, show that the anisotropy in the local interaction (A) does not change the 
qualitative behavior of the necklace model on the cubic lattices with D=l-3. In order to study the effect of a nonzero 
anisotropy S ^ in the 'itinerant' part (~ t) in Eq. Q it is therefore sufficient to consider only the isotropic case of 
the local Kondo interactions, i.e., J x = J z = J, (A = 1). 
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FIG. 7: The effect of <5 on the behavior of gap versus t for the three dimensional cubic lattice. The gap vanishes at the critical 
point with different scaling behavior for < 8 < 1 and 8 > 1, (8 = 2, 3, 5). 



When 6 is nonzero, we find three excitation modes in Ea. l|ll|l which now all show dispersion. The minimum of 
excitations which introduces the energy gap appears again at the antiferromagnetic reciprocal wave vector Q. fn 
contrast to the previous case, the excitation in z-direction (w z ) is important and defines the gap for the easy-axis type 
anisotropy 8 > 1. The minimum excitation energies for the modes at Q are given by 



For < 8 < 1, we have Eg < Eg. However, at 8 = 1 the energy gap of all excitations coincide. If we consider 5 > 1 
the situation is reversed and Eg < Eg. This is important because the vanishing of gap and the appearance of soft 
modes define the transition from Kondo singlet to the antiferromagnetic phase. 

In the case of one dimensional space (D = 1) the energy gap is always nonzero for < 8 < 1. The dependence 
of the gap on the anisotropy of the 'itinerant' term (5) in Eq.Q is smooth as shown in in Fig|SJ However, there is 
an abrupt change in the value of gap at 8 = 1 which is related to the contribution of the third soft mode (w z ) to 
the lowest energy part of the model. The situation is almost the same for 8 > 1 where the gap never vanishes and 
therefore no phase transition appears for increasing t. 

The effect of 8 anisotropy in the two dimensional lattice is plotted in Fig|5{a). As discussed earlier the gap E g is 
defined by the Eg and E v g for < S < 1 and Eg for i5 > 1, (8 = 2, 3, 5). For 8 < 1 the effect of anisotropy is weak and 
it changes the critical point (t c ) slightly. The dependence of the critical point (t c ) is not a simple function as will be 
given in the next section. The qualitative behavior is the same for all 6 < 1 and the gap exponent is approximately 
v ~ 1. We observe a jump in the critical point which appears at 8 = 1. Because at this point another soft mode 
(lu z ) is added to the excitations at the critical point which needs a larger hopping strength t to overcome them. If 
we describe the gap close to the critical point by E g = A(8)\t — t c \ u the coefficient A(S = 1) is almost half the value 
for the other cases (A(8 ^ 1)). When 8 > l,t c is reduced and the Kondo singlet phase is limited to a smaller region. 
In this case the interaction between the itinerant electron (tS) dominates the Kondo term (J) for smaller values of t. 
The change of universality in approaching the critical point is also obvious from FigUJa) . 



J, / . ts . 
-D( -), 



(18) 



10 



Increasing t beyond t c leads to antiferromagnetic order as discussed in Sect.E] For comparison we have also plotted 
the AF order parameter (staggered magnetization) versus t in Fig|SJb). The staggered magnetization is zero in the 
singlet phase until t = t c . For high values of t (t 3> J), the antiferromagnetic order parameter is close to saturation. 
The reduction of t increases the effect of the local exchange interaction which favors the singlet phase. At the critical 
point [t — t c ) the staggered magnetization is destroyed by the AF fluctuations connected with the soft modes. This 
happens exactly at the point where the Kondo singlet energy gap mentioned in previous paragraph vanishes. 

We also show the effect of nonzero 8 on the three dimensional cubic lattice in Fig[7] As discussed earlier, the 
minimum of excitations in x- or y-direction gives the scale of energy (Eg , E^) for < i5 < 1 and for 5 > 1. The 
effect of 8 on the behavior of energy gap versus t is weaker than the case of D = 2. The gap vanishes at the critical 
point which has a slight dependence on 5. There is a jump on the critical point at 8 = 1 which is the result of three 
soft modes at this point while for < 8 < 1 there are two soft modes and for 8 > 1 only one soft mode. Although 
8 is a global term in the interaction there is no change in the scaling behavior of the gap close to the critical point. 
The gap exponent for < 8 < 1 is the same as for zero 8, namely v ~ 0.55. However, we have observed different 
behavior for 8 > 1 where the gap vanishes close to critical point more rapidly. But we were not able to find a scaling 
exponent in this regime because here t c varies too strongly with 8. In the other cases we have determined the scaling 
exponents numerically from the selfconsistent solutions and found them to lie close to v = i for 2D and v = 1 for 3D. 
The deviations are mostly due to the numerical inaccuracy of t c determination. This is the expected result for the 
present mean field approach. In the special case (8, A) = (0,1) we have expanded the selfconsistent equations Eq. l|14l) 
around E g =0 which indeed leads to the analytical mean field scaling behaviour for E g . 



V. ANTIFERROMAGNETIC PHASE 



The mean field approach can be simply extended to the antiferromagnetic phase for D = 2,3. For large enough tj J 
the long range antiferromagnetic order sets up which means the condensation of a component of local spin triplets. We 
assume the condensation in the x-component of the spin triplet (tk, x ) at the antiferromagnetic wave vector (Q = (it, 7r) 
for D=2), 



tk,x = V Nt8 k ^ Q + T]k, x ■ (19) 

In Ea. Q19[l t is the mean value of the x-component of the spin triplet in the ground state and rjk^ is its quantum 
fluctuations. We replace t^ lX into the Fourier trasformed of Eq.© to derive the mean field Hamiltonian assuming 
(sfc) = s and (tk,x) — i. The mean field Hamiltonian will be in the following form, 

H AF = E$ F + X>z( fc )*L*M + XXWfafcV* + iljk, y ), (20) 

fc k 

where O represent the new operators after Bogoliubov transformation. The ground state energy (Eq F ) is defined 
below, 

E$ F = E + NP(n + ^-Dts 2 ), (21) 

where Eq has been defined in Ea. (|12fl . The ground state energy Eq F is minimized with respect to /i, s and t which 
leads to the following equations, 



•s 



■ 2 5 J x + J z 1 v J 7 (fc) 1 v X% 
4 4Dt 2N ^ V D AN ^ ujJk) 

k k 



i2 _ 5 Jx + Jz 1 \ ^ 1 1 \ ^ -X"- <22) 

~n2-^ I 77TT AN 2^1 ,,,JhV ^ ' 



Furthermore, 



4 ADt 2N ^ + 4A ^ co z (k) ' 



X% = I(J x -^ )( l + ^) + ^s 2 (l + ^), 

Xt = ^(J x -J z )(l-S^-)+Dts 2 . (23) 
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FIG. 8: Dependence of the quantum critical point t c (in units of 3 X ) for D — 2 on itinerant (S, full line) and local (1-A, dashed 
line) anisotropy parameters. The lines without symbols are from the closed analytical form in Eg. 12511 . the curves with symbols 
are obtained from numerical solution Eg. 1221 or Eg. 1141 . 



Generally the equation for s should be solved selfconsistently and then /j, and i can be found directly from the above 
equations. However, for A = 1, the above equation (Ea. (|22[l ) for s is simplified to a direct integral expression. If in 
addition 5 = the result in RefQ is recovered. For the bipartite lattice which we have assumed the antifcrromagnetic 
order parameter is defined on A- or B-sublattices as, 

(S X ) A = -{S X ) B = si, (24) 

Here m s =si is the staggered moment. It is plotted for D = 2 and different 5 values in Fig.|jjb). The appearance of 
long range antiferromagnetic order (nonzero m s ) defines the quantum critical point t = t c . This is identical to the 
value where the Kondo singlet gap vanishes as seen in the conjugate plot in Fig.|SJa). The formulation presented in 
this section allows us to obtain an explicit expression for the critical point (t c ) for the special cases (S, A) — (0, A) and 
(5, A)= (S, l)which is given by 

Here we used again j Xj y(k) = j(k) and J z (k) = 6~f(k). For the other cases one has to obtain t c numerically by the 
selfconsistent solution of Ea. H22l) or equivalently by Ea. (|14fl . For D=2 we have plotted the value of t c versus 1 — A 
(at fixed (5 = 0) or versus S (at fixed A = 1) in Fig. [S] The bare lines correspond to the analytical solution of Ea. i|25|) 
while the lines with symbols correspond to the numerical solution for t c . One notices that on approaching (5, A) — > 
(1,1) where i c =0.88 (the two uppermost points in Fig.|SJ| the slope of t c (A) or t c (S) diverges. This signifies the change 
of universality class from U(l) to SU(2) at the point (S, A)=(l,l) which leads to the apperance of a third soft mode 



VI. SUMMARY AND DISCUSSIONS 



We have studied the quantum critical properties of the fully anisotropic Kondo-necklacc model in dimension D=l-3. 
We used the bond-operator representation of spin variables with constraints implemented in the mean field approxi- 
mation. For general anisotropy ratios (5, A) we have calculated the critical values (t/ J) c and the scaling exponenets 
of the Kondo singlet gap close to the quantum critical point which correspond to the mean field exponents. 

We find that in general the anisotropy A of the local Kondo term influences the qualitative behavior of the singlet 
gap only little, although the value of the critical ratio (t/J) c depends strongly on A. If A decreases the Kondo singlet 
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is destabilized and the critical value (t/J) c decreases. On the other hand the anisotropy of the inter-site term which 
mimics the itinerant electrons has a weak effect on the position of the critical point (Fig. |SJ and leaves the universal 
behavior similar as in the conventional XY-type (S = 0) Kondo necklace model. On approaching the SU(2) Heisenberg 
case (5, A) = (l,l) t c exhibits a singular behaviour as function of anisotropics. This is most pronounced in D=2 where 
also the scaling coefficient (not the exponent) of the gap close to the QCP changes considerably as compared to the 
5 = case. Furthermore the value of the critical (t/J) c with properly defined scale J exhibits a jump at 5=1. This 
peculiar effect is due to a mode crossing of the x,y and z branches of excitations at the AF wave vector as function 
of 6. 

We have also derived and solved the selfconsistency equations on the magnetic side for general (5, A). This allows 
us to give an explicit expression for the quantum critical point value t c as function of anisotropy parameters (5, A) 
for the cases 6=0 or A=l for dimension D=2,3. To discuss the scaling exponents beyond mean field approximation 
the renormalisation of triplet excitation energies close to the QCP caused by fluctuations in the bosonic variables has 
to be included. 
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VIII. APPENDIX 



The explicit form of Hi and H% introduced in Sect. ITTll in terms of bond operators are as follows: 

H 2 = ((*3i, 1/ *n,«-/i.c.)(*^ lV W _ft - a ) + (4,xW-^0(*kxW _ ' l - c -) 

(n,m) 

+ £(4, a *n,3/ ~ h.c.){tl n x t m , y - h.c.)j . (26) 

The terms appearing in Hi describe the interaction between triplets. Since the triplet occupation is very small in the 
ground state, the effect of Hi on the mean field result is negligible. 

The remaining part of H t is expressed by H% which gives zero contribution to the ground state energy in the mean 
field approximation. 



H U 
H 3 - - 

(n,m) 



n,m) 

+ S^S^tn^z + ^n !Z s n)(^L,a:^"i,a — ^ln,y^rn,x)] + h.C^j . (27) 
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